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Preparing a quantum system in a pure state is ultimately limited by the nature of the system's evolution in 
the presence of its environment and by the initial state of the environment itself. We show that, when the 
system and environment are initially uncorrected and arbitrary joint unitary dynamics is allowed, the 
system may be purified up to a certain (possibly arbitrarily small) threshold if and only if its environment, 
either natural or engineered, contains a "virtual subsystem" which has the same dimension and is in a state 
with the desired purity. Beside providing a unified understanding of quantum purification dynamics in 
terms of a "generalized swap process," our results shed light on the significance of a no-go theorem for exact 
ground- state cooling, as well as on the quantum resources needed for achieving an intended purification 
task. 

Cooling of quantum systems toward their ground state plays a central role across low-temperature physics 
and quantum science, by providing the key to unlock novel phases of matter and quantum behavior - as 
exemplified in settings ranging from laser cooling of atoms and molecules 1 " 3 to dynamical nuclear polar- 
ization in solid- and liquid-state nuclear magnetic resonance 4 ' 5 , and cooling of mechanical resonators 6 " 10 . From a 
quantum control standpoint, the task of cooling (or "refrigeration," in the language of quantum thermodyn- 
amics 11 ) may be viewed as an instance of dissipative pure- state preparation, which is in turn closely related to the 
more general task of purification - namely, the ability to steer the system from an arbitrary initial state to a final 
state with higher purity. Within quantum information processing (QIP), access to pure states is presumed in all 
quantum computation models that can provably achieve an exponential speed-up over classical ones 1213 , and cold 
ancilla qubits are critical to the success of fault-tolerant quantum error correction 14 . As a result, schemes for 
cooling and purification are being actively investigated 15 " 17 , and underlying assumptions and implications for- 
malized with added rigor 18 " 21 . 

While in practice a variety of system-dependent imperfections and technological constraints will inevitably 
hinder the achievable performance, a fundamental question is to determine what ultimate limitations may 
nevertheless exist on the sole basis of some generic, system-independent assumptions on the underlying dynamics. 
Specifically, assume that arbitrary unitary evolution is allowed on the target system S together with its envir- 
onment E, starting from arbitrary factorized initial conditions. To what extent does the initial, typically highly- 
mixed state of £, limit the degree of purity attainable on S in principle 7 . Conversely, if the environment E and its 
initial state can be controllably engineered, what are the minimal resources for purification (cooling) of S to be 
guaranteed to a prescribed accuracy? 

Our main contribution in this work is the identification of necessary and sufficient conditions for exact as well as 
approximate purification and ground-state cooling, given the above ideal scenario. Our starting point is a trivial 
example: if both S and E are two-dimensional systems (qubits), purification of S is clearly possible in principle 
only if the initial state of E has a lower entropy, in which case the optimal purification dynamics simply amounts 
to swapping the two initial states. In a general open -system setting, our strategy is to make precise the intuition 
that purity can still only be exchanged but not created between subsystems, albeit the latter need no longer 
coincide with the natural ones. The relevant notion is provided by the concept of a "virtual" subsystem as & factor 
of a subspace of a larger state space, as introduced by Knill et al. 22 in the context of quantum error correction and 
extensively used in QIP 23 " 28 . 

Our results complement existing work and advance current understanding in several ways. While a no-go 
theorem for ground- state cooling under initial system-thermal bath factorization was recently established in 
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Figure 1 | The system of interest, S, may be generally coupled to a 
quantum bath, 5, and an engineered auxiliary system, A. We collectively 
refer to the pair (B, A) as the environment. The initial state on 
TLse = Hs®He = Hs® (Hb ® Ha) is assumed to be fully factorized with 
respect to this partition, i.e., p S E = Ps ® Pe = Ps ® (Ab ® Pa) - The joint 
dynamics is generated by a total Hamiltonian of the form 
H = H 0 + H c (t) = (Hs®Ie+1s®H e + H S e) +H c (f), where the control 
Hamiltonian H c (t)= ue(t)H c> £ acts trivially on B. If dim(HsE) < °o> 
complete propagator controllability is ensured in the generic case where 
the Lie algebra of skew- symmetric operators generated by the control 
Hamiltonians { iH c j], together with the natural "drift" Hamiltonian iH 0 , is 
the whole $u(ds x ds) 35 . If so, there exist some time T > 0 and control 
functions Uf(t), t e. [0, T], that allow to reach any element in H(Hse) to 
arbitrary precision. For our discussion, it is not essential to specify how the 
control actions are enacted. For instance, if Ha — C, our setting includes 
open-loop control of S via a semiclassical controller 36 ' 37 . In this case, B is 
controlled via its interaction with S, yet indirect controllability of B given an 
arbitrary initial state of S still suffices for complete joint controllability, as 
we assume 38 . If dim (7-^) > 1> dynamics in the presence of a coherent 
"quantum controller" and/or an an engineered quantum reservoir 16 ' 37 may 
be accounted for. In this case, the uncontrollable component B may couple 
to both S and A in general. 

Ref. 20, our analysis further clarifies that such a no-go strictly applies 
only to exact cooling. Notwithstanding initial factorization, no fun- 
damental limit exists to arbitrarily accurate purification and ground- 
state cooling in principle, so long as the environment is effectively 
infinite- dimensional, and capable of supporting a sufficiently pure 
virtual subsystem. From a quantum-simulation standpoint, this rein- 
force the conclusion that a simulated ancillary environment consist- 
ing of a single qubit suffices for enacting arbitrary open-system 
dynamics, so long as it can be measured and reset to a sufficiently 
pure state 29 , as recently demonstrated in trapped-ion experi- 
ments 30,31 . Conceptually, our analysis points to a generalized swap 
process as the unifying physical mechanism through which any puri- 
fication or groundstate cooling dynamics may ensue from joint unit- 
ary evolution, as opposed to known special instances limited to small 
dimension and/or a fixed (thermal) initial environment state 20 ' 21 . 
From an open- system quantum- control perspective, our general pic- 
ture may be exploited to design procedures for purification and 
ground-state cooling via environment (or "reservoir") engineering, 
as potentially relevant to a growing number of quantum technolo- 
gies, see e.g. Refs. 29,32,33 and references therein. Interestingly, 
within quantum foundations, our results have also implications for 
dynamical reduction models 34 : in order for the "wave-function col- 
lapse" predicted by the standard von Neumann postulates to be 
consistently reproduced by underlying open-system dynamics, the 
environment interacting with the system must, again, harbor a suffi- 
ciently pure virtual subsystem. 

Results 

Setting. The general setting we consider is depicted in Fig. 1. The 
target quantum system S, with associated Hilbert space Hs of 
dimension d s , is coupled to a quantum environment £, with 
associated Hilbert space He of dimension d E , which may generally 
include both a component that is not directly controllable (a physical 
"bath", B) and a fully controllable auxiliary system (or "ancilla", A). 



We take d E > d s , so that we may decompose d E = d s d F + d R , with d F 
being is the integer part of d E /d s , and d R < d s the rest. While we 
further assume that d s < °° in what follows, we may formally extend 
our results to infinite- dimensional target systems of interest (notably, 
quantum oscillators) by imposing a finite-energy constraint. 

A key assumption is that no correlations are initially present 
between the constituents of the joint system, i.e., the initial state is 
factorized, p SE = Ps ® Pe> with p E a trace- class operator in case 
d E = oo. Other than that, and unlike in Refs. 20,21, no restriction is 
placed on either p s or p E which, in particular, need not be thermal. 
We shall denote by U ; (p £ )} the eigenvalues of p E , considered with 
their multiplicity and in non-increasing order. 

While the inclusion of both a bath and an ancillary system allows 
for different physical scenarios to be discussed within the same 
framework (see caption), the central mathematical assumption is 
that suitable Hamiltonian control is available on S + E together, so 
that any unitary operator in H(Hse) can be obtained at some time T. 
In control-theoretic terms, this is equivalent to assuming complete 
joint propagator controllability 35,39 . Hence, at any given time T, the 
joint evolution of S + E is described by some U SE (T) that we are free 
to choose. The conditions for this to be possible have been extensively 
investigated within the geometric control framework. At least if Hse 
has finite dimension, complete controllability is generic 35 , albeit effi- 
cient constructive methods for control design are still object of 
ongoing research, along with controllability conditions for infinite- 
dimensional quantum systems 40,41 . 

Starting from factorized initial conditions, the reduced state of the 
system after the unitary (controlled) evolution takes place is given by 

p's = Tr £ (p'se) = Tr £ ( U SE p s ® Pe Use) • ( 1 ) 

Exact purification of S is attained if p' s is pure irrespective of the 
initial state p s , that is, p' s =\\l/)(\//\ for some \\j/)eHs, so that 
Tr(ps 2 ) = 1. However, this requirement is too strong in practical 
situations of interest. We say that (s-)approximate purification of S 
can be attained if the state of S may be brought to within distance s 
from a pure state irrespectively of the initial p s , that is, there exists 
\\l/)eH s such that 

d(p' s ,WM)<8, Vp s . (2) 

Here, d(X,Y)=^Tr(\X-Y\)='^\\X-Y\\ 1 is the quantum total- 
variation distance, which is a natural measure of distinguishability 
between quantum states 12,26,27 . (Note that the robustness requirement 
with respect to the system initialization makes our purification 
notion stronger than used in both Refs. 20,21. For given initial p s 
and p E with known spectrum, an upper bound on the purity of the 
final state p's maybe additionally established 21 .) Exact purification is 
recovered by requesting s = 0. In the following, we shall consider 

0<8<1. 

The fact that the joint dynamics Pse^Pse is unitary is equivalent 
to the preservation of the spectrum of the joint density operator at 
any time. However, one still intuitively expects purification of a 
"portion" of the overall system to be possible in an appropriate sense, 
the limitations on what can be achieved stemming from the initial 
state of E. Let us first consider a trivial example. 

Example 1- Suppose that both the target system and the envir- 
onment are a qubit. The factorized initial state can then be parame- 
trized by the maximum eigenvalue of its two components p s , p E , say, 
1/2 < p s , p E < 1 respectively, with the value 1/2 corresponding to a 
fully mixed state. That is, p SE = diag(p s , 1 — p s ) (x) diag(p £ , 1 — p E ). 
Since for qubits the von Neumann's entropy S{p) is completely 
determined by, and is a decreasing function of, the maximum eigen- 
value of the state, we can pursue a direct information-theoretic ana- 
lysis. Achieving maximal purification is thus equivalent to achieving 
the (reduced) state p's in Eq. (1) with minimum entropy with respect 
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to the choice of U S e- Using the standard definitions of joint and 
conditional entropy 12 , we may write 

S(p' s )=S{S)=S(S,E)-S{E\S), 

where S(S,E) =S(p'se)=S(Ps®Pe)> ana " the conditional entropy is 
maximal when the state is factorized. Hence, the maximal purifica- 
tion is attained by either swapping the states (when p E > p s ), or 
leaving them as they are (when p E < p s ). In other words, some 
purification is possible only if the entropy of the auxiliary qubit is 
lesser than the one of the system qubit, and exact purification can 
only be achieved if the former is in a pure state to begin with. 

Despite its simplicity, this example suggests a general strategy to 
tackle the purification problem: given a target system to be purified, if 
in its environment we may identify a "subsystem" of the same 
dimension, that is initially in a more pure state, all we need to do is 
to swap these two states. Formalizing this idea leads to the rigorous 
conditions we are seeking. 

Main result: necessary and sufficient conditions for purification. 

In common physical situations, subsystems may be naturally 
identified with (distinguishable) quantum particles and/or degrees 
of freedom, and their state space directly associated to different 
factors of the overall tensor-product Hilbert space. This view is 
not, however, sufficiently general to capture all relevant settings 
that arise both physically and in the context of QIP applications. 
Within quantum error correction theory, for instance, "noise- 
protected" quantum-information-carrying logical degrees of free- 
dom are associated with virtual subsystems that typically do not 
correspond with the original qubit subsystems 22 ' 25 ' 42 . This more 
general subsystem notion will also be key to our analysis. 
Mathematically, a virtual quantum subsystem S of a larger system 
E (the environment in our case) is associated with a tensor factor H~ s 
of a subspace of He 22 ' 25 , that is, we may write 

n E = {H~ s ®H F )eH R , (3) 

for some factor H F and a (generally non -trivial) remainder space Hr. 
System E is said to be initialized in subsystem S if its state may be 
decomposed as p E = p~ s ®p F @0 R , where 0 R is the zero operator on 
Hr and p F a state on H F \ in particular, E is initialized in a subsystem 
pure state if p~ s = \(p)((p\, for \g?)eH~ s 26 ' 43 ' 44 . While virtual subsystems 
are most compactly described in terms of an operator-algebraic 
characterization 22 ' 23,27,28 , a basis with the correct tensor/direct 
product structure may also be straightforwardly constructed (see 
Methods). We are now ready to state our central result: 

Theorem. Assume complete unitary controllability and factorized 
initial conditions p s ® Pe and on Hs®H F - Then the following 
conditions hold: 

(i) For every s > 0, s-approximate purification of S may be 
achieved if there exists a decomposition of He as in Eq. (3), with 
H~ s ^Hs> and a pure-state initialization of E in S, 
Pe = \<p)(<p\®Pf® sucn tnat 

d(p E ,p E )<£. (4) 

(ii) Exact purification ofS(e = 0) may be achieved if and only if 
the initial state of the environment has exactly the form 
Pe = \<P) (<p\® Pf®Qr> for some \(p)eH~ s . 

(iii) ^-approximate purification is always possible provided that 
s>s(p E ), where 

d F 

~s( Pe ) = ~s = 1-J2¥Pe)>0- (5) 
;'=i 

z-purification is optimal whenever d R = 0. In particular, arbit- 
rarily accurate purification ( 8 = 0, 8 > 0) is always possible for 

d E = oo 



Part (i) of the above theorem can be easily proven by considering a 
unitary operator W SE that at some time T swaps the state of S with the 
one of its isomorphic copy S, which is initially in a pure state \qj) {q>\. 
With the precise definition of W SE being given in the Methods sec- 
tion, the basic observation is to note that if p E satisfies Eq. (4), then it 
can be written as 

p E = p E + Ap E , ±Tr(|Ap £ |) < £ . (6) 
By implementing the swap dynamics, it thus follows that 

p' s = Tr E (w SE p s ®p E wl E ^ 

= Tr B [|^>(^|®(p s ®p J7 eo fi ) + W S BPs®Ap B wy ( 7 ) 

= \f)((p\+8(p s ®kp E ), 

where 8 is a trace-preserving completely-positive map and hence a 
trace-norm contraction 12 . Since, together with Eq. (6), this implies 
that 

d{p's,\(p){(p\)<z, Vp s , 

the target system S is g-purified, as desired. It is immediate to see that 
the same argument also applies if s = 0. In other words, the condition 
of Eq. (4) is always sufficient for s -purification with s > 0, indepen- 
dently of the dimension and the initial state of E. 

Establishing that Eq. (4) remains necessary is relatively straight- 
forward for exact purification [as in part (ii)], but more subtle in the 
approximate case [part (iii)]. While full proofs are given in the 
Methods section, the gist of the argument showing why ^purifica- 
tion is indeed always possible for 8 > e may be summarized as follows. 
Assume that for an initial state p SE = p s 0 Pe> the desired purifica- 
tion can be attained at some final time T. Then there exists an 
orthogonal projector, say, II r = \\l/)(\//\s ® I E , such that 
Tr(n T ps £ ) > 1 — g, for all p s . If we define a new projector 
n 0 = uI e T1 t Use> this condition clearly also implies that Tr(n 0 ps£) 
> 1 — g. This inequality shows that a pure subsystem of dimension d s 
may be identified to within distance s from the initial jo int state as 
well. The tricky part is to establish that this in turn implies the 
existence of an g-pure subsystem in the environment alone. 

In order to do so, we may consider the worst-case scenario, that is, 
a fully mixed (infinite temperature) initial state on S, with 
p SE = p SE = (l/ds)Is®p E - The idea is to construct a projector of 
the form n 0 = /<j®ni, where Hi is a projector on d F eigenvectors 
of p E with highest eigenvalues, which projects on a subspace, say 
Hi £ He, of the same dimension of IT r . This is the best possible 
strategy whenever d R = 0, and we may show that: 

Tr(n lP£ )=Tr(n 0 p S£ ) = l-g, ~8<8. (8) 

Accordingly, the subspace Hi, onto which II x projects, collects 
(1 — e) of the total probability. The existence of such a subspace 
may be shown to be equivalent to the existence of a virtual subsystem 
S, such that E is 8-close to pure-state initialization in S, as desired. 

Our theorem points to an interesting dichotomy between finite - 
vs. infinite- dimensional environments. If d E < oo, g -purification of S 
may or may not be achievable, depending on whether the conditions 
on the spectrum of p E imposed by Eq. (8) are fulfilled, for arbitrary 
p s . If d E = oo, however, then 8 = 0 and for any trace-class state of E 
and any fixed s > 0, a sufficiently pure subsystem always exists. We 
illustrate how to explicitly construct such a £-pure subsystem in the 
case where the target system is a qubit, as the generalization to a 
higher- dimensional system (qudit) is straightforward. 

Let p E be a trace-class environment state, and consider its spectral 
representation, say, p E = x pj \j) (j\ > 5^/- iPj = ^' ^ e i( ^ en " 
tification of the desired £-pure subsystem may be accomplished by 
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Figure 2 | The target system (with <i s -dimensional state space Hs) is 
coupled to an infinite- dimensional quantum bath (with state space Hb), 
initially in an arbitrary state p B . To construct a subsystem of B which is 
arbitrarily (yet not perfectly) pure, we identify a finite-dimensional 
subspace H\ that collects the first M eigenvectors of p B accounting for ( 1 — 
e) of the total probability. To complete this virtual subsystem, we only need 
to identify (d s — 1) orthogonal subspaces H, each of dimension M. 
Purification is then attained by swapping the virtual subsystem's state with 
the one of the target system. 

identifying two orthogonal subspaces Hi, H2 £ He each of dimen- 
sion M, one of which accounts for (at least) (1 — s) probability. Since 
p E is trace class, hence its spectrum is absolutely summable, for any 
e > 0 there exists an M large enough such that >M Pj < 8 - Define 
H\ =span{|j)} ;=1 M , and H2 any M-dimensional subspace 
orthogonal to Hi. From these two subspaces, we can easily con- 
struct a subsystem decomposition as in Eq. (3), with (H~ s )=2, 
dirndl?) = M, such that the final reduced state p's is £-close to a pure 
state. The strategy is pictorially illustrated in Fig. 2. The general qudit 
case can be obtained along the same lines, by considering d s copies of 
the M-dimensional subspace, where again only one accounts for (at 
least) (1 — e) of the total probability. Similar considerations also 
apply to typical physical scenarios where d F ~d E /ds~dE^>ds, in 
which case nearly arbitrary accuracy £—0 may still be achieved in 
principle. 

Our results show how there is no fundamental limit to arbitrarily 
accurate purification when coupling the target system to an effec- 
tively infinite-dimensional environment. Exact purification, on the 
other hand, would require a sufficiently large number of eigenvalues 
of p E to be precisely zero. Since this is not a generic condition, in 
particular it cannot be obeyed if p E is thermal, the no-go theorem of 
Ref. 20 is recovered. With this general conceptual framework in 
hand, we next proceed to examine in more detail relevant applica- 
tions, beginning from the special important case of ground-state 
cooling. 

Ground-state cooling given initial system-bath factorization. 

Consider a setting where, as in Fig. 2, the environment consists of 
a physical bath (E = B), and let H s denote the (free) Hamiltonian of 
the target system S, so that the corresponding initial energy is 
Tr(H sPs ). 

Assume first that the minimum eigenvalue £ min of H s is not degen- 
erate, in which case exact cooling of S to its ground state entails 
preparing it in the unique pure state |i/f gs ) corresponding to eigen- 
value £ min . It is then a straightforward corollary of our theorem that 
exact ground-state cooling can be obtained only if the environment 
contains a virtual subsystem of the same dimension of the target, 
which is initialized in a pure state. Under the complete joint unitary 



controllability assumption, however, the ability to prepare a given 
pure state also imply the ability to prepare any pure state in Hs. 
Hence, the existence of a pure virtual subsytem of the environment 
is also necessary for exact cooling, fully consistent with the conclu- 
sions reached in Ref. 20. 

On the other hand, suppose that only e-approximate purification 
may be achieved in the sense of Eq. (2), so that the state of S can only 
be cooled down to within distance s > 0 from the unique ground 
state \4f gs ) of H s . Then the final energy of the system may be esti- 
mated as 

Tr [H s p' s ] = Tr [ff s (( 1 - e) | ^ gs ) | + 8T ex ) ] 

— (1 ^)^min H~ £^max? 

where T ex and £ max denote some state in the orthogonal complement 
to the ground manifold and the maximal eigenvalue of H s , respect- 
ively. Accordingly, approximate ground-state cooling may be 
attained with an "excess" energy that is upper-bounded by sE max . 
We already observed that when E is infinite-dimensional, £ can in 
principle be chosen arbitrarily small, albeit not zero. Thus, as soon as 
one allows for approximate yet arbitrarily good cooling, the no-go 
theorem can be effectively evaded 20 . 

If £min has degeneracy d gs > 1, being able to prepare a pure state 
still suffices for exact ground- state cooling, but is no longer needed. 
Sufficient and necessary conditions for approximate cooling in a 
degenerate subspace may be derived using the same reasoning used 
in establishing necessity of our condition for s > s - by finding a 
virtual d s -dimensional subsystem S such that E is g-close to initiali- 
zation in a subspace of dimension d gs in H~ s . 

Arbitrary purification and ground-state cooling with an engineered 
qubit reservoir. From an open-system quantum-control perspective, 
our theorem may be used to explicitly characterize what quantum 
resources may suffice to arbitrarily purify/cool the target system, by 
coupling it to a suitably engineered environment (E = A). Let us focus 
on the simplest yet paradigmatic case in which S is a single qubit and A 

consists itself of N qubits, so that Ha — (C 2 ) x . 

Building on the previous discussion, identifying the desired virtual 
qubit- subsystem entails to split Ha into two isomorphic, orthogonal 
subspaces. For the resulting "virtual state" to be approximately pure, 
we further require the probability for A to be found in one of such 
subspaces to be much higher than the one for its complementary. A 
natural approach is to invoke a "typical subspace" argument. Let 
each auxiliary qubit be prepared in the same state, say, p = dmg(q, 
1 — q), 1/2 < q < 1, with respect to a standard basis {|0), 1 1)}, so that 
the joint initial state p SA = p s ® p® N . As AT grows, the state of A will 
populate with increasing probability the e-typical subspace. Recall 
that a sequence x(N) of N zeroes and ones, in which each entry is 
chosen independently at random with probability P(0)=q, 
P(l) = l-q, is e-typical if 12 

2 -N(S(x) + e) < p < 2 -N(S(x)-e) ? 

or, equivalently, its total Shannon entropy is e- close to N times the 
binary entropy of the single symbol. Let T(N,e) be the set of e-typical 
sequences. In the quantum case, such a set naturally generalizes to 
the e-typical subspace: in our qubit setting, the latter is spanned by 
those computational basis states that include (approximately) qN 
zeroes and (1 — q) N ones: 

H r{ N,e)^^n{\x(N))\x(N)eT(N,e)}. 

Let now Uj-(N,e) denote the orthogonal projection onto the e-typ- 
ical subspace. Then the following asymptotic result holds (see e.g. 
Theorem 6.3 in Ref. 45): 
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Jim qtJp (N) = lim Tr(n T(N>e) p® N ) = 1. (9) 

Furthermore, for any fixed e > 0 and a sufficiently large N, the size 
of the typical subspace satisfies: 

dim(H T{N ,e))<2 N(S(p)+e) . 

Hence, if e is sufficiently small, the dimension of the e-typical 
subspace becomes less or equal than half of the total space dimension 
as soon as NS(p) <N— 1, or, «S(p) < 1 — l/N. Therefore, provided 
that the entropy of each of the auxiliary qubits is strictly less than one, 

namely p # - 1, the typical subspace's dimension will become less 

than half of the dimension of Ha for large enough N. If so, we know 
how to explicitly construct a unitary transformation W S a that 
achieves (optimal) £ -purification in principle: it suffices to swap 
the state of S with the state of a virtual qubit system S that exploits 
the typical- subspace structure. We further illustrate this strategy by 
specializing, again, to ground-state cooling. 

Example 2- Assume, similar to Example 1, that the initial state of 
the target system p s = diag(p s , 1 — p s ), with respect to the qubit 
energy basis, say, = {|*Ag S )> l*Aex)} andp s < q. The action of the 

unitary transformation W S a may be explicitly described by intro- 
ducing a factorized basis ® \j(N))} on TLs®Ha> where 
{[/'(N)> « [/typ), i/ntyp)} in the large-iV limit, with {|j typ )} and {\j ntyp )} 
denoting orthonormal bases for the typical subspace and its ortho- 
gonal complement, respectively. The idea is then to swap ~2 NS ^ 
typical basis states which have non-zero probability and are assoc- 
iated to |^ex)> with ~2 N<s ( p ) non-typical basis states which are in 
tensor product to | i^ gs ) but are associated to low probability. If we 
compute the final energy of the system, by using Eq. (9) we obtain 
Tr[Hsps]~gtyp(N)£ m in = (1 — s)£min> with arbitrarily small e (hence 
e) as N — > oo, as desired. 

Altogether, our results imply that, for a target qubit system, arbit- 
rary accuracy in purification and cooling may be achieved through 
fully coherent (unitary) interaction with sufficiently many copies of 
any auxiliary qubit state which is not the completely mixed one. It is 
interesting to notice, however, that repeated interactions with an 
identically prepared qubit do not suffice in general: the generalized 
swap operation needs to simultaneously operate on multiple qubits 
of the engineered environment, pointing to an intrinsic non- 
Markovian action. 

Robust pure-state preparation with finite control iterations. As a 

final application of our framework, our main theorem may be used to 
understand and characterize the control resources involved in a 
stronger form of purification, whereby the goal is to bring the state 
of S to a predetermined target pure state |*A) target 6 ^' not necessar ily 
related to the system's ground state - so-called "global asymptotic 
stabilization" in control-theoretic parlance 33 ' 35 ' 39 ' 44 . In particular, the 
case where | »A)target is an entangled pure state on a multipartite n- 
qubit target system provides an important quantum-stabilization 
benchmark. While it is well-known that access to a single resettable 
ancillary qubit A, along with complete unitary control over S and 
fully coherent "conditional" interactions between A and S, suffices to 
engineer arbitrary dynamics on S 29 ' 37 (and hence achieve the desired 
stabilization task) in principle, our result sheds light on the 
thermodynamical foundation of this result. With reference to the 
general setting of Fig. 1, suppose for simplicity that no uncon- 
trollable bath is coupled to S (H SB = 0), and that B represents the 
physical degrees of freedom which enact, possibly together with 
coherent control between S and A, the resetting process on A. 
Then, in order for stabilization of S to be achievable with arbitrary 
accuracy 8 starting from an arbitrary environment state p E = 



Pb®Pa> an effectively infinite-dimensional environment is 
necessary. Furthermore, exact pure-state stabilization is only 
achievable provided that A may be perfectly refreshed, which in 
turn requires B to be perfectly initialized in a pure, two-level 
virtual subsystem. Remarkably, if these conditions are met, an 
arbitrary n-qubit pure state | iA)tar ge t ma y in fact be dissipatively 
prepared by using a finite number, n, of suitably defined control 
iterations 46 . 

Experimentally, controlled dissipation mechanisms are becoming 
available in a growing number of scalable platforms for universal 
"digital" open-system quantum simulators, including trapped- 
ion 30,31 and superconducting qubit technologies 47 . In the above- 
mentioned experiments on 40 Ca + ions, for instance, the required 
re-initialization dynamics of the ancilla qubit to a reference pure 
state was realized through a combination of coherent control on A, 
in conjunction with optical pumping followed by spontaneous emis- 
sion. While a number of details are important and require careful 
consideration in practice, conceptually it is this step that ultimately 
grants access to virtual subsystems whose states are sufficiently pure, 
and can thus be swapped with those of the physical degrees of free- 
dom to be purified and/or cooled. 

Discussion 

We have identified sufficient conditions for purification and ground- 
state cooling of a quantum system of interest to be achievable in 
principle, under the two assumptions of initial system- environment 
factorization and complete unitary controllability. Such conditions 
are also necessary in most realistic situations, where the environment 
is much larger than the target system. While in essence these condi- 
tions make rigorous an intuition that is both compelling and natural 
in retrospect - namely, that purity can only be "swapped" across 
appropriately defined quantum subsystems - we have shown how 
these conditions allow to both elucidate fundamental limitations in 
harnessing open -system dynamics as well as identify new opportun- 
ities for control engineering. In particular, our analysis makes it clear 
that arbitrarily accurate purification and/or ground- state cooling is 
always possible in principle as long as the relevant environment is 
effectively infinite- dimensional, with a no-go result 20 only emerging 
in the limiting case of zero error. 

From a control-theory standpoint, an interesting direction for 
further study is to characterize what (more stringent) limitations 
on quantum purification and cooling may arise upon relaxing the 
assumption of complete controllability for S + E. We envision that 
the existence of a sufficiently pure virtual subsystem in the envir- 
onment will still be a necessary and sufficient condition, albeit iden- 
tification of the relevant subsystem structure will be carried out in 
this case by exploiting the dynamical- symmetry decomposition 
associated to the reachable control sub-algebra, in analogy to dynam- 
ical error-control strategies and encoded tensoriality in QIP 48 49 . 

Lastly, it is interesting to comment on our results in relationship to 
the third law of thermodynamics in its dynamical formulation - the 
so-called "unattainability principle", namely, the impossibility to 
cool a system to absolute zero temperature in finite time 11 . 
Throughout our discussion, we have deliberately made no explicit 
statement on the time T needed to implement the required general- 
ized swap transformation W SE (T). For a standard thermodynamic 
setting where the bath is given, and is initially in a generic trace- class 
state (say, thermal at non-zero temperature), we have showed that 
arbitrarily small cooling error, s > 0, may be achieved only if a 
sufficiently large subspace of the bath can correspondingly account 
for less than s probability. This, in turn, translates into an increas- 
ingly complex (energetically "delocalized") action of the swap trans- 
formation W S e(T) to be implemented. Since realistic control 
Hamiltonians are inevitably constrained (e.g., bounded in amplitude 
and/or speed, as stressed in Refs. 15,16), the limit of perfect accuracy, 
s — > 0, can only be approached asymptotically in time, T— » oo. While 
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this supports the validity of the third law under typical conditions, it 
is our hope that our general subsystem-based approach may prove 
useful to deepen our understanding of fundamental performance 
bounds in more complex thermodynamic scenarios, including 
"quantum-enhanced" refrigeration as recently proposed in Ref. 50. 

Methods 

Subsystem construction and generalized swap operation. Starting from a general d- 
dimensional Hilbert space H, a "virtual subsystem structure" as used in the main text 
can be identified by constructing a basis with the correct tensor/direct sum structure. 
The main steps may be summarized as follows: 

• Identify a (d x X d 2 ) -dimensional subspace Hi, 2 , so that H^Hi, 2 ®Hr, where 

Hr = Hi 2 =HQHi ; 2- 

• Inside such a subspace, choose d x mutually-orthogonal subspaces H h 2, each of 
dimension d 2 , so that we may decompose Hi,2 — ®f =x H h 2- 

• Pick an orthonormal basis in each of the summands, say, { | <ft) k , k= 1 , . . . , d 2 }, for 
j = 1, d\. We can then establish the following identification: 

|^), = l^>s®l^>^ H lf2 ^H s ®H F , 

and obtain the desired subsystem structure, with dim(TCs) = d\, dim(H F ) = d 2 , 
respectively. 

Consider now, specifically, a subsystem structure as given in Eq. (3) on the 
environment Hilbert space, namely, H F = (H~ s ®H f )®Hr, and let {|^ ; >s}> { |<Ws}> 
{|<^)f}> {|XmM be orthonormal (ordered) bases for Hs, H~ s , H F , Hr, respectively. We 
may define the required generalized swap unitary operator W SF through its action on 
the element of an orthonormal basis. That is, consider the (ordered) basis of Hs®H F 
given by: 

{\^) s ®\h)s®\^)F}^{\^) s ®\x m ) R }, 

for all j, k, €, m. The action of Wse is then defined by: 

v^B(l^> s e|z«> Jl ) = l* / > s e|zm>i l - 

Proof the main theorem. Assume that, as in the main text, we write d E = d s d F + d Ry 
with d R < d s ^ d E , and let p SF — Ps ® Pe denote an arbitrary joint initial state on 
H S ®H E . 

Proof of part (ii). The fact that the existence of an e-pure subsystem in the envir- 
onment suffices for e-purification (g > 0) has already been proved in the text. We 
show here that for the case of exact purification (e = 0), the existence of a purely- 
initialized, <i s -dimensional subsystem in He is indeed also necessary. 

Recall that exact purification is equivalent to the existence of an orthogonal pro- 
jector, II r = \\\j){^j\ s ® I E , such that Tr(nrp' S£ ) = 1, and that upon defining 
n 0 = uI e IItUse, this also implies that 

Tr(n 0 p SE )=Tr(n 0 p s (8)p E ) = l, Vp s . (10) 

This in particular means that the support of p s ® Pe is included in the range of Il 0 . Let 
us consider d s specially chosen initial states p s , associated to an orthonormal basis 
[\4>j>£ ofH s ,thatis, 

{Ps®P E )j = Wj)($j\s®PE, j=h--- A- 

Since their supports are mutually orthogonal, and each of them has dimension 
rank(p £ ), it follows that: 

rank(n 0 )> ^rank(|^)(4|(g)p £ ) =d s rank(p £ ). (11) 

;=1 

On the other hand, since rank(Il r ) = d E , we also have rank(n 0 ) = d E . Together with 
Eq. (11), this implies: 

rank(p £ )<^, 
d s 

and hence, being an integer, rank(p £ ) < d F . Call Hi = supp(p E ) <= He, and construct a 
d s -dimensional virtual subsystem of He as described above. By construction, p E is 
purely initialized in the first elements of the basis associated to the d s -dimensional 
subsystem Hs, leading to the desired conclusion. 



Proof of part (Hi). Let us define the following two quantities [see also Eq. (5)]: 

d F 

~s(p E )=~s=l- J2^(Pe)> 
j=i 

^R(p E )= £ R=^r^d F +i(p E )- 

d s 

We next proceed to show that: 

1 . A lower bound e 0 exists for purification of S, independently of the initial state p s ; 

2. Purification up to e = e 0 + s R is always possible by properly identifying a sub- 
system in He alone and then swapping it with the target. 

1. Determining e 0 — We look for necessary conditions on e > 0, so that s-purifica- 
tion of S can be attained at time T by some joint unitary transformation U SF - 
Again, this means that there exists an orthogonal projector, II r = I^X'Als ® h, 
such that Tr(n r p' S£ ) > 1 — £, for all p s . Upon defining U 0 = uI e IItUse as 
above, this also implies that 

Tr(n 0 p S£ ) = Tr(n 0 p s (x)p £ )>l-£, Vp s . (12) 

Thus, a pure subsystem of dimension dg may be identified to within e-distance 
from the joint initial state as well. While Eq. (12) must hold for all p$, in order to 
determine the desired lower bound we consider a worst- case scenario where p s 
= (l/d s )Is and, correspondingly, the initial joint state p SE = p SE = (l/ds)ls®p E - 

In fact, consider a basis in which p SE is diagonal, ordered in such a way that its 
eigenvalues are non-increasing. The eigenvalues of p SE are the eigenvalues of p E , 
each multiplied by (l/d s ) and repeated d s times. Given that n 0 has rank d E , the 
maximal purification achievable in this case correspond to n 0 projecting on the 
first d E eigenvalues. It is then easy to show that: 

J2h(Ps E ) = y2^j(pE) + -r^d F +i(pE)- ( 13 ) 

Since, to guarantee e-purification, the (i £ -ranked projector n 0 must satisfy Eq. 
(12) in particular for p SE = p SE , Eq. (13) implies the following lower bound e 0 : 

£>£ 0 =£-£r. (14) 

We remark that so far nothing guarantees that £ 0 -purification is attainable for 
any initial state. 

2. Affa/m'ftg£-purification.- From Eq. (13), we infer that there exists a subspace H\ 
of He alone, with dimension d F , that accounts for 1 — £ = 1 — £ 0 — 8r of the trace 
of p E . We can thus consider the subspace H\ ^ H E that collects only the one- 
dimensional eigenspaces corresponding to the first d F eigenvectors of p E . The 
last step is to start from Hi to construct a virtual subsystem S, such that E is £- 
close to pure- state initialization in S. 

We can in fact identify additional (d s — 1) orthogonal subspaces in He, say, 
{l~Lj}j S =2 , all isomorphic to Hi and composed of eigenspaces of p E , so that, by 
following the general subsystem construction described above, have: 

He = (jl Hj^j ®H R ^ (H~ s ®H f ) ®H r , 

where H s ^ H s , d\m(H F ) = d Fy and H R = H E G ®f = iHj, dim(H R ) = d R . Let n x be 
the orthogonal projector onto Hi, and define n 0 =/s®IIi. By construction, n 0 has 
rank d s d F < d E . It thus follows that: 

Tr(n 0 p S£ ) = rr(n lP£ ) = i-£. 

Now notice that with respect to the subsystem decomposition above, we may write 
111 = \(p)((p\®I F ®0 R for some \(p)eH~ s , and 

Pe = Pe + a Pe = \v) (<?l ®*f ®0 r + Ap E , 

with t f = - — rdiag(/li (p E ), ■ ■ ■ Ad F (Pe))- Accordingly, with respect to the decom- 
position He = Hi ®H±, we may write Ap E = Ap l ®Ap^, with 

Api = Y^diag(^i(p £ ), • • • M(Pe))> 

Apj 1 =diag(A di7 + i(p E ), . . . M e (Pe))- 

Since these matrices correspond to the positive and negative- semidefinite part of 
Ap £ , it follows that 
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l -Tr(\Ap E \)= l - [-Tr(A Pl )+Tr(Ap^)] =~s. 



We may thus conclude that p E admits a e-pure subsystem, and by using the gen- 
eralized swapping we can guarantee ^-purification of the target, as claimed. 

Note that whenever d R = 0, we have s = s 0 and thus our generalized swap operator 
attains the best possible purification. If, additionally, d E = °°, this also formally 
corresponds to d F = °° hence e = 0. We have then explicitly shown in the main text 
how to achieve purification up to arbitrary finite accuracy s > 0. 
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